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In this talk I review studies of hadron structure functions in bosonizcd chiral quark 
models. Such models require regularization and I show that the two-fold Pauli— 
Villars regularization scheme not only fully regularizes the effective action but 
also leads the scaling laws for structure functions. This scheme is consistent with 
other computations of the pion structure function in that model. For the nucleon 
structure functions the present approach serves to determine the regularization 
prescription for structure functions whose leading moments are not given by matrix 
elements of local operators. Some numerical results are presented for the spin 
structure functions and the role of strange quarks is addressed. 



1. Introduction 

In this talk I review the computation of hadron structure functions in the 
Nambu-Jona-Lasino (NJL) model 1 . This is a particularly simple model for 
quark interactions with the important feature that the quarks can be inte- 
grated out in favor of meson fields 2 . The resulting effective action for these 
mesons possesses soliton solutions 3 . According to the large- Nq picture 4 
of Quantum-Chromo-Dynamics (QCD) these solutions are interpreted as 
baryons. Quantization of the soliton then also yields baryon wave-functions 
in such meson models. The construction of hadron wave-functions is not 
possible in QCD. This represents the main obstacle for the computation of 
hadron properties from first principles. As the model adopts the symmetry 
properties of QCD, the current operators in the model correspond to those 
of QCD. As a consequence matrix elements of the current operators as com- 
puted in the model are physical. In particular it is interesting to analyze 
the hadronic tensor that parameterizes the deep-inelastic-scattering (DIS) 
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and compare the model predictions with empirical data. This picture has 
led to interesting studies of hadron structure functions in bosonized chiral 
quark models. In this talk I will present the results of refs. 5,6 ' 7,8 . These 
studies build up a consistent approach by computing the hadronic tensor 
(or equivalently the forward virtual Compton amplitude) from the gauged 
meson action. For the nucleon structure functions similar studies have 
been reported in refs. 9 ' 10 ' 11 . There no attempt to compute the structure 
functions from the gauged action was made but rather the less convincing 
assumption was made that the model predictions for the constituent quark 
distribution can be identified with QCD quark distributions. However, I 
refer to those articles for a more expatiated presentation of numerical re- 
sults. 

This talk is organized as follows. In Section 2 I introduce the NJL 
model as an effective meson theory and utilize pion properties to determine 
the model parameters. Section 3 describes the subtleties for extracting the 
structure functions that arise in this model from regularization. The pion 
structure function is considered as an example. In Section 4 I review the 
construction of baryon states in the soliton picture. The following Section 
sketches the computation of nucleon matrix elements of the hadronic tensor 
and the extraction of the structure functions in the Bjorken limit. Finally 
in Section 6 I present some numerical results for the spin structure func- 
tions gi and gi and compare them to experimental data by means of the 
transformation to the infinite momentum frame and subsequent DGLAP 
evolution. Section 7 serves as a short summary. 

2. The NJL Model for Chiral Dynamics 

The NJL model is a quark model with a chirally invariant four quark inter- 
action. Semiclassical bosonization is achieved via introduction of effective 
meson fields for the quark bilinears in that interaction. Then the quark 
fields can be integrated out by functional methods. This yields an effective 
action for meson degrees of freedom, 



A[S, P] = -iN c Tt K \og [i$ - (S + h 5 P)\ - — / d 4 xtrV(S, P) . (1) 



Here V is a local potential respectively for the effective scalar and pseu- 
doscalar fields S and P that are matrices in flavor space. In the NJL model 
it reads V = S 2 + P 2 + 2rfioS with too being the current quark mass matrix. 
Since the interaction is mediated by flavor degrees of freedom, the number 
of colors is merely a multiplicative quantity. The functional trace (Tr) orig- 
inates from integrating out the quarks and induces a non-local interaction 
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for S and P. For simplicity I will only consider the isospin limit for up (u) 
and down (d) quarks: mo. u = fno,d = m a- 

A major concern in regularizing the functional (1), as indicated by the 
cut-off A, is to maintain the chiral anomaly. This is achieved by splitting 
this functional into 75-even and odd pieces and only regulate the former, 

N 2 

Tr A log [i? -(S + i l5 P)] = J] c 4 Tr log [— DD 5 + A 2 - ie] 

i=0 

-z^Trlog[-D(D 5 )- 1 -ze] , (2) 
with iD = i$ - (S + ij 5 P) and i~D 5 = -i@ - (S - i<y 5 P) . (3) 

The double Pauli-Villars regularization renders the functional (1) finite 
with co = 1, Ao = 0, Y^ii=o c i = 0- The 75 _ odd piece is conditionally 
finite and not regularizing it properly reproduces the chiral anomaly. For 
sufficiently large G one obtains the VEV, (S) = ml that parameterizes the 
dynamical chiral symmetry breaking, from the gap-equation, 



2 _ ( m _ mo ) = -AiN c m z2 Ci J J^j* ^ + m " + A ? ~ ie \ ■ ( 4 ) 



Substituting S — (S) — ml in cq. (1) shows that to plays the role of a mass 
and is referred to as the constituent quark mass. 

In the next step I utilize pion properties to fix the model parameters 
and introduce the pion field ir via 

S + iPj 5 = m (Uy< 5 = mcxp (*^75 7? • t) . (5) 

Sandwiching the axial current between the vacuum and a single pion state 
yields the pion decay constant f n — 93MeV in terms of the polarization 
function II (q 2 , x), 



f 

Jo 



L = 4:N c mg / dxU(m n ,x) 



d 4 k 



U(q 2 ,x) = J2 Cl 7^a 1 [ fc2 + a; ( 1 _ • T ) m - -m 2 - A 2 + ie] , (6) 
i=o ^ n ' 1 

where to t = 138MeV is the pion mass. The Yukawa coupling constant, g, is 
determined by the requirement that the pion propagator has unit residuum, 

7 = mc ^l! dx[mlu{mlx)] - (7) 
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In the chiral limit (m„ = 0) this simplifies to f„ = m/g. Finally the current 
quark mass is fixed from the condition that the pole of the pion propagator 
is exactly at the pion mass, 

m =4:NcmGml / dx U(ml , x) . (8) 
Jo 

It is also worthwhile to mention that expanding eqs. (2) and (5) to linear 
and quadratic order in ir and w M , respectively yields the proper result for 
the anomalous decay 7r° — > 77. This is the direct consequence of not 
regularizing the 75-odd piece. 

Before discussing nucleons as solitons of the bosonized action (1) and 
the respective structure functions it will be illuminating to first consider 
DIS off pions. 



3. The Compton Tensor and Pion Structure Function 

DIS off hadrons is parameterized by the hadronic tensor W^(q) where q 
is the momentum transmitted from the photon to the hadron. 

The tensor W^ v (q) is obtained from the hadron matrix clement of the 
current commutator by Fourier transformation and is parameterized in 
terms of form factors that multiply the allowed Lorentz structures. By per- 
tinent projection of the hadron tensor these form factors can be extracted. 
Finally the structure functions are the leading twist contributions of the 
form factors. These contributions are obtained from computing (q) in 
the Bjorkcn limit: q 2 — ► —00 with x = —q 2 /p ■ q fixed. 

An essential feature of bosonized quark models is that the derivative 
term in (1) is formally identical to that of a non-interacting quark model. 
Hence the current operator is given as J M = qQ'f^q, with Q a flavor ma- 
trix. Expectation values of currents are computed by introducing pertinent 
sources in eq (2) 

iD — > iH + Qi> and iD 5 — > iT> 5 - Qj> (9) 

and taking appropriate derivatives. In bosonized quark models it is con- 
venient to start from the absorptive part of the forward virtual Compton 
amplitude* 1 

T^(q) = J d^e^<{p, S \T(J^Or(0))\p,.s} 
W^(q) = ±-\m[T^(q)}, (10) 

a The momentum of the hadron is called p and its spin eventually s. 
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because the time-ordered product is straightforwardly obtained from 

T (JW(0)) = 5v J* 5vM Tr A log & - (S + , 75 P) + Qf\ , (ID 

as defined from eq. (2) with the substitution (9). 

Pion-DIS is characterized by a single structure function, F(x). For 
its computation the pion matrix element in the Compton amplitude (10) 
must be specified. For virtual pion-photon scattering it is obtained by 
expanding cqs. (2) and (5) to second order in both, tt and v^. Due to 
the separation into D and D5 this calculation differs considerably from the 
simple evaluation of the 'handbag' diagram. For example, isospin violating 
and dimension-five operators appear for the action (2). Fortunately all 
isospin violating pieces cancel yielding 

F{x) = \{AN c ?)^[m%U{ml,xj\ , < x < 1 . (12) 

The same result is obtained by formal treatment of the divergent handbag 
diagram and ad hoc regularization 12 . The cancellation of the isospin vio- 
lating pieces is a feature of the Bjorken limit: insertions of the pion field 
on the propagator carrying the infinitely large photon momentum can be 
safely ignored. Furthermore this propagator can be taken to be the one for 
non-interacting massless fermions. This implies that also the Pauli-Villars 
cut-offs can be omitted for this propagator. That, in turn, leads to the 
desired scaling behavior of the structure function in this model and is a 
particular feature of the Pauli-Villars regularization. A priori it is not ob- 
vious for an arbitrary regularization scheme that terms of the form q ' / Af 
drop out in the Bjorken limit. 

From eqs (7) and (12) it is obvious that F(x) = 5/9 for < x < 1 in 
the chiral limit (m, =0). It must be noted that this refers to the structure 
function at the (low) energy scale of the model. To compare with empirical 
data, that are at a higher energy scale, the DGLAG program of perturbative 
QCD has to be applied to F(x). Such studies show good agreement with 
the data 13 . 



4. The Nucleon as a Chiral Soliton 

Solitons are a non-perturbative and stationary configurations of the meson 
fields. To determine that configuration for the meson theory (1) I consider 
the hedgehog ansatz 

Un{r) = exp(ir • fF(r)) and (C/h(^)) 75 = exp (i^f ■ fF(r)) 
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for the pion field (5). The corresponding single particle Dirac Hamiltonian 
reads 

h = a ■ p + f3 m [cosF + ij5T ■ f sini* 1 ] . (13) 

Evaluating the action functional (1) in the cigcnbasis of h gives the energy 
functional 14 in terms of the eigenvalues, e a , 

E[F] = ^ (1 - sign( ev) ) ev - ^ E E c 4 " V^+^j 

+m 2 Jl j d 3 r (1-cosF) (14) 



for a baryon number one configuration. Here V denotes the unique quark 
level that is strongly bound by the soliton and are the eigenvalues 
when the soliton is absent. The soliton profile F(r) is then obtained from 
cxtrcmizing E self-consistently 3 . 

States possessing good spin and isospin quantum numbers are generated 
by taking the zero-modes to be time dependent 15 

U(f,t)=A(t)U H (f)A\t), (15) 

which introduces the collective coordinates A(t) € SU{2). The action func- 
tional is expanded 16 up to quadratic order in the angular velocities 

if-ti = 2A i (t)A(t). (16) 

The coefficient of the quadratic 13 term defines the moment of inertia , a 2 [F]. 
Nucleon states \N) are obtained by canonical quantization of the collective 
coordinates, A(t). This is analogous to quantizing a rigid rotator. The 
eigenfunctions are the well-known Wigner D-functions 

with 73 and J3 being respectively the isospin and spin projection quantum 
numbers of the nucleon. This then allows to compute matrix elements of 
operators in the space of the collective coordinates 15 : 

(N\^ti(T a A^T b A)\N) = -l(N\I a J b \N} and ti = -J/a 2 [F]. (18) 



b A liner term does not arise due to isospin symmetry. 

c Functional integrals are evaluated using the eigenfunctions <p a of the Dirac Hamilto- 
nian (13) in the background of the chiral angle F(r). Thus all quantities - like the 
moment of inertia - turn into functionals of F(r). 
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For later use I note that the valence quark wave-function receives a first 
order perturbation 

* v (r , t) = e-^Ait) { Mr) + \ £ M?) W ' m \ ■ (19) 

The moment of inertia, a 2 [F] is order Nc, thus, upon quantization (18), 
this rotational correction is subleading in 1/Nc- 

Also, later in this talk I will generalize the soliton approach to the case 
of three flavors. The quantization of the soliton proceeds along the same 
line by taking A(t) e SU(3). However, there is one essential difference: 
ST/ (3) is only an approximate symmetry because the current mass of the 
strange quark differs substantially from that of the up and down quarks. 
This will be reflected in the Hamiltonian for the collective coordinates A(t) 
by an explicit flavor symmetry breaking contribution. Nevertheless this 
Hamiltonian can still be diagonalized exactly by numerical means 17 . As 
a result the eigenfunctions are not simple Wigner D-functions of specific 
S£/(3) representations but rather linear combinations thereof. These linear 
combinations reflect the fact that spin-i baryons are no longer pure octet 
states but acquire admixture of higher dimensional SU (3) representations. 
For example the nucleon becomes 

\N) = | A, 8) + Cl [F] | A, 10) + c 2 [F] | A, 27) + . . . 

m s + m m s + m 

where the constants of proportionality, Cj are computed from the self- 
consistent soliton. The exact wave-functions will be employed when study- 
ing nucleon properties in the three flavor model. For further details on 
quantizing the soliton in flavor SU(3) I refer to the reviews 18 . 

5. Nucleon Structure Functions 

DIS off nucleons is described by four structure functions: Fi(x) and F 2 {x) 
are insensitive to the nucleon spin. Those associated with the components 
of the hadronic tensor that contain the nucleon spin are g\(x) and g 2 (x). 

As argued in section 3, the quark propagator with the infinite photon 
momentum should be taken to be the one for free and massless fermions. 
Thus, it is sufficient to differentiate (Here D and D 5 are those of eq (3), 
i.e. with = 0.) 

^ E c * Tr { (- DDs + A *) _1 [ q2 ^ ^y 1 ^ D ^ ^y 1 ^ Q2 } } 

+ ^Tr {(-DDg)- 1 [Q 2 i, (Vy 1 i,T> 5 + Dtf ^y 1 ^) 5 Q 2 ] } , (20) 
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with respect to the photon field v^. I have introduced the (. . .)5 description 

to account for the unconventional appearance of axial sources in D5, cf. 
ref. 8 . Substituting (15) for the meson fields that are contained in D and D 5 , 
computing the functional trace up to subleading order in 1/Nq using a basis 
of quark states obtained from the Dirac Hamiltonian (13), yields analytical 
results for the structure functions. I refer to ref. 8 for detailed formulae for 
other structure functions and the verification of the sum rules that relate 
integrals over the structure functions to static nucleon properties. As an 
example I restrain myself to list the contribution to g\ (x) which is leading 
order in 1/Nc- 

d\ 
2tt' 



itoX 



\i=0 



+^(0T3(l-a 3 ) 7 5^(C-Ae 3 )e^ A , (21) 



where the subscript (P) indicates the pole term. 

Before discussing numerical results I would like to mention the unex- 
pected result that the structure function entering the Gottfried sum rule is 
related to the 75-odd piece of the action and hence does not undergo reg- 
ularization. This is surprising because in the parton model this structure 
function differs from the one associated with the Adler sum rule only by the 
sign of the anti-quark distribution. The latter structure function, however, 
gets regularized in the present model, in agreement with the quantization 
rules for the collective coordinates. 



6. Numerical Results for Spin Structure Functions 

Unfortunately numerical results for the full structure functions in the double 
Pauli-Villars rcgularization scheme, i.e. including the properly regularized 
vacuum piece are not yet available. However, in the Pauli Villars rcgu- 
larization the axial charges are saturated to 95% or more by their valence 
quark (19) contributions once the self-consistent soliton is substituted. This 
provides sufficient justification to consider the valence quark contribution 
to the polarized structure functions as a reliable approximation since e.g. 
the zeroth moment of the leading structure function gi is nothing but the 
axial current matrix element. This valence quark level is that of the chiral 
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soliton model and its contributions to the structure functions should not 
be confused with valence quark distributions in parton models. In general, 
it should be stressed that the present model calculation yields structure 
functions, i.e. quantities that parameterize the hadronic tensor, but not 
(anti)-quark distributions. The latter would require the identification of 
model degrees of freedom with those in QCD. However, here only the sym- 
metries (namely the chiral symmetry) and thus the current operators in the 
hadronic tensor are identified. 

As in the case for the pion, the model calculation yields the nucleon 
structure function at a low energy scale. In addition it must be noted 
that the soliton is a localized object and thus not a representation of the 
Poincare group. As a result the computation of structure functions is not 
frame-independent. It is appropriate to choose the infinite momentum 
frame (IMF) not only because it makes contact with the parton model but 
also because it is that frame in which the support of the structure functions 
is limited to the physical regime < x < 1. Choosing the IMF amounts to 
the transformation 19,20 



where /rf(^) denotes the structure function as computed in the nucleon 
rest frame. So the program is two-stage, first the transformation of the 
model structure function to the IMF according to eq (22) and subsequently 
the DGLAP evolution program 21 to incorporate the lnjg 2 corrections. In 
the current context it is appropriate to restrain oneself to the leading order 
(in a s ) in the evolution program because higher orders require to iden- 
tify quark and antiquark distributions in the parton models sense. In the 
present model calculation this is not possible without further assumptions 01 . 
The low energy scale, Qq = 0.4GeV 2 , at which the model is assumed to 
approximate QCD has been estimated in ref. 5 from a best fit to the ex- 
perimental data of the unpolarized structure function, -£2(2;). The same 
boundary value is taken to evolve the model prediction for polarized struc- 
ture function, gi(x), in the IMF (22) to the scale Q 2 of several GeV 2 at 
which the experimental data are available. For the structure function g 2 (x) 
the situation is a bit more complicated. First the twist-2 piece must be 
separated according to 22 



Jimf(x) 



1 



/ RF (-ln(l-x)), 



(22) 



l-x 




(23) 



We assume, however, that the gluon distribution is zero at the model scale. 
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Figure 1. Model predictions for the polarized proton structure functions xg\ (left panel) and 
x 92 (right panel). The curves labeled 'RF' denote the results as obtained from the valence 
quark contribution to (20). These undergo a projection to the infinite momentum frame 
'IMF' (22) and a leading order 'LO' DGLAP evolution 21 . Data are from SLAC-E143 24 . 

and evolved analogously to gi(x) (which also is twist-2). The remainder, 
92 (x) — g^ w (x) is twist-3 and is evolved according to the large-iVc scheme 
of rcf. 23 . Finally, the two pieces are again put together at the end-point 
of the evolution, Q 2 . In figure 1 I compare the model predictions for the 
linearly independent polarized structure functions of the proton to exper- 
imental data . In figure 2 I compare the model predictions for both the 
proton and the neutron (in form of the deutcron) not only to the recently 
accumulated data but also to other model predictions. Surprisingly the 
twist-2 truncation, i.e. eq (23) with the data for gi(x) at the right hand 
side gives the most accurate description of the data. However, also the 
chiral soliton model predictions reproduce the data well. Bag model pre- 
dictions have a less pronounced structure. 

As mentioned in section 4, the chiral soliton model can be generalized to 
three flavors. Appropriate projection gives a prediction for the strangeness 
contribution to structure functions 7 . In figure 3 this contribution is shown 
for the spin structure function g\ (x) . The interesting feature is that gf^ (x) 
has both positive and negative pieces. This is a nice example showing that 
a small As = Jq dx g[ s ^ (x) (strange quark contribution to the nuclcon spin) 
does not imply that strangeness structure function itself is small. 

7. Conclusions 

I have presented studies of the nucleon spin structure functions in chiral 
soliton models. For this purpose I considered the bosonized NJL model as 
a simplified model for the quark flavor dynamics. Although the bosonized 
version is an effective meson theory, it has the interesting feature that the 
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Figure 2. Model predictions for the polarized proton structure functions xg^ for proton 
and neutron (deuteron) and comparison with data from E143 24 (open diamond) and E155 25 
(open square) and their combination (solid circle). The full line is the twist-2 truncation (23) 
of data for gi (x). Dashed-dotted 26 and dotted 27 lines are bag model calculations, the short 
dashed lines represent the present chiral soliton model 6 and long dashed line that of ref. 11 . 
(This is a slightly modified figure from ref. 25 .) 




Figure 3. The strange quark contribution to the nucleon structure function gi(x) in 
the infinite momentum frame (IMF) and in leading order (LO) evolution from the model 
scale Ql = 0.4GcV 2 to Q 2 = 3GcV 2 . 



quark degrees of freedom can be traced. This is very helpful for considering 
structure functions. It turned out that additional correlations are intro- 
duced due to the unavoidable regularization which is imposed in a way to 
respect the chiral anomaly. Hence a consistent extraction of the nucleon 
structure functions from the Compton amplitude in the Bjorken limit leads 
to expressions that are quite different from those obtained by an ad hoc 
regularization of quark distributions in the same model. I also showed that 
within a reliable approximation the numerical results for the spin dependent 
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structure functions agree reasonably well with the empirical data. 

I would like to thank the organizers for this worthwhile workshop. The con- 
tributions of my colleagues L. Gamberg, H. Reinhardt, E. Ruiz Arriola and O. 
Schroder to this work are gratefully acknowledged. This work has been supported 
by the Deutsche Forschungsgemeinschaft under contract We 1254/3-2. 
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